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ABSTRACT 

In this article we present the results of the fourth order perturb ation a nalysi s of the 
gra vitational metri c theo ry of gravity /(%) = x 3 ^ 2 ? developed by iBernal et al. ( 2011b ) 

and iMendoza et al.l (|2013l) . The theory accounts in detail for the mass from the observations 
of 12 Chandra X-ray clusters of galaxies, without the need of dark matter. The dynamical ob¬ 
servations can be obtained in terms of the metric coefficients of the metric theory of gravity up 
to the fourth order of approximation, in perturbations of v/c. In this sense, we calculate the first 
relativistic correction of the theory, which is relevant at the outer regions of clusters of galaxies, 
in order to reproduce the observations. 


Subject headings: gravitation - relativity - galaxy clusters 


1. Introduction 

The observations of type la supernovae, the 
anisotropies observed in the microwave back¬ 
ground, the acoustic oscillations in the bary- 
onic matter, the power-law spectrum of galaxies 
and gravitational lenses among others, represent 
strong evidences for the standard cosmological 
model, the so-called ACDM concordance model 
From recent observations of the European space 


mission Planck , the contribution of the baryonic 
matter to the present content of the matter-energy 
density of the universe was inferred to be only 
5%, while the dark sector constitutes ~ 95%, of 
which 27% is Cold Dark Matter (CDM) and 68% 
is dark e nergy or a positive cosmolo gical constant 
A ( Planck Collaboration et aOl2013l) . 


The dark matter component was postulated in 
order to explain the observed rotation curves of 
spiral galaxies, as well as the mass to light ratios 
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in giant galaxies and clusters of galaxies, the ob¬ 
served gravitational lenses and the structure for¬ 
mation in the early universe, among other astro- 
physical and cosmological phenomena. On the 
other hand, the dark energy or a cosmological con¬ 
stant has been postulated to ex plain the acceler¬ 
ated expansion of the universe (IPerlmutter et al. 

199§- 

The ACDM model adjusts quite well most of 
these observations. However, direct or indirect 
search of dark matter candidates has yielded null 
results. In addition, the lack of any further evi¬ 
dence for dark energy opens up the possibility that 
there are no dark entities in the universe but in¬ 
stead, the theory associated to these astrophysical 
and cosmological phenomena needs to be modified. 

Current models of dark matter and dark energy 
are based on the assumption that Newtonian grav¬ 
ity and Einstein’s general relativity are valid at all 
scales. However, their validity has only been ver¬ 
ified with high precision for systems which scales 
are no larger than the Solar System one. In that 
sense, is conceivable that both, the accelerated ex¬ 
pansion of the universe and the stronger gravita¬ 
tional force required in different systems, represent 
a change in our understanding of the gravitational 
interactions. 


From the geometrical point of view, modi¬ 
fied theories of gravity are viable alternatives to 
solve the astrophysical and cosmological prob¬ 
lems that dark matter and dark ener gy are try¬ 


ing to solve (see e. g. Schimming fc Schmidt 2004 


Nojl 

n & Odintsov 1 

2011 

). In this sense, 


2011: ICapozziello &; Faraonil 


ity which attempts to supplant the dark compo¬ 
nents of the Universe, must account for two crucial 
observations: the dynamics observed for massive 
particles and the observations of the deflection 
of light for mass l ess pa rticles. As extensively de¬ 
scribed by IWilll ( 1993 ). when working with the 
weak field limit of a relativistic theory of gravity 
in a static spherically symmetric spacetime, the 
dynamics of massive particles determine the func¬ 
tional form of the time component of the metric, 
while the deflection of light d eterm ines the form of 
the radial one (see also Willi 2006 . and references 
therein). 


The first successful modification in the non- 
relativistic regime to deal with these issues was 
the Modified Newtonian Dynamics (MOND) (see 


Famaev fc McGaughl 2012 . for a review). Due 
to its phenomenological nature and its success in 
the non-relativistic limit, it is understood that 
any fundamental theory of modified gravity should 
adapt to it on galactic scales in the low accel¬ 
erations regime. However, from the study of 
groups and clusters of galaxies it has been shown 
that, even in the deep MOND regime, a dominant 
dark matter component is required in these sys¬ 
tems^ ( 60 t o 8 0% of the dynamical or virial mass). 
Angus et al.l (12008 ) showed that the central region 
of galaxy clusters could be explained with a halo 
of neutrinos with mass of 2 eV (which is about the 
value of the experimental upper limit). But on the 
scale of groups of galaxies, the central contribution 
cannot be explained by a contribution of neutri¬ 
nos with that mass. Moreover, MOND/AQUAL, 
the Lagrangian formulation of MOND, is not able 
to reproduce the observed gr avitation al l ensing 
for different system s (see e.g. iTakahashi fe Chiba 


2007 : iNataraian fc Zhaql2008[ ). mainly because it 
is a non-relativistic description, and as such it can¬ 
not explain gravitational lensing and cosmological 
phenomena, which require a relativistic theory of 
gravity. 

Through the years, there have been some at¬ 
tempts to find the relativistic extension of MOND. 
The fi rst su ccessf ul attempt was proposed by 
Bekensteir] ( 2004 1 who formulated a Tensor- 
Vector-Scalar (TeVeS) theory. This approach 
presents some cumbersome mathematical com¬ 
plications and it cannot rep roduce crucial astro - 
physical phenomena (see e.g. lFerreras et al.ll2009h. 
Another recent construction bv lDemir & Karahan 
(12014 propose a relativistic version of MOND 
through modifications in the energy-momentum 
tensor. In an empirical way, they recover the 
MONDian limit modifying the dynamics sector. 

As will be discussed in se ction |3l the extended 
gravity theory proposed by iBernal et al.l ( 2011bh 
is equivalent to the MONDian description in some 
systems, for example in spherical s ymmetric ones 


but w ith remarkable advantages. Men doza et al 


( 2013 ) have performed a second order perturba¬ 
tion analysis of the gravitational me tric theory of 
gravity /(%) = y 3 / 2 proposed first by Bernal et al 
( 2011 b). They have shown that the theory ac¬ 


counts in detail for both observational facts: it is 
possible to recover the phenomenology of flat rota¬ 
tion curves and the baryonic Tully-Fisher relation 
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(i.e. a MOND i an-like weak field limit; se e e.g. 
Mikroml Il983at iFamaev fc McGaughl I2012L and 
references therein), which includes the mass of the 
system in the gravitational field’s action; and the 
details of observations of gravitational lensing in 
individual, groups and clusters of galaxies, with¬ 
out the need of any dark matter component. 

In the present work, we extend the perturba¬ 
tion analysis up to the fourth order of the the¬ 
ory, in powers of v/c, where v is the velocity 
of the sys tem and c i s the speed of light. As 
shown by Sadeh et al.l (120151) and first hypothe¬ 
sized bv fWoitak et al.l ( 201 li b there exist observa¬ 
tional relativistic effects of the velocity of galax¬ 
ies in the edge of galaxy clusters, showing a dif¬ 
ference of the inferred background potential with 
the galaxy inferred potential. With this motiva¬ 
tion in mind, we have calculated fourth order rel¬ 
ativistic corrections t o the f( y) = y 3 / 2 theory 
( Bernal et al . 2011b : IMendoza et al.l 112013 ). and 
show that they can account to fit the observations 
of 12 Chandra X-ray clus ters of galaxies presented 
bv IVikhlinin et all (2 006 ). 

The article is organized as follows. In section [2] 
the weak field limit for a static spherically sym¬ 
metric metric of any theory of gravity is estab¬ 
lished and we define the orders of perturbation 
to be used throughout the article. In section [3] 
we show the particul ar metric theory proposed by 
Bernal et al. ( 2011b ) to be tested with the astro¬ 
nomical observations. The results from the per¬ 
turbation theory for the vacuum field equations 
up to the fourth order in perturbation for such 
metric theory are presented. With these results, 
we obtain the gravitational acceleration generated 
by a point-mass source and its generalization for 
extended systems, particularly for applications to 
clusters of galaxies. In section 0] we establish the 
way to fit the corrections to the metric coefficients 
from the observations of 12 Chandra X-ray clus¬ 
ters of galaxies. Finally, in section [5] we show the 
parameters calibration of the generalized gravita¬ 
tional acceleration obtained in section [3] through 
the observationally derived dynamical mass of the 
clusters of galaxies, and we discuss our results. 


2. Perturbations in spherical symmetry 


to the metr ic theory developed b y iBernal et al 


( 2011b ) and IMendoza et all ( 20131 ). Many of the 
results developed in t his section were obtained by 


Mendoza et al.l ( 2013 ). This section is kept here 


for completeness of the article. 


Perturbations applied to metric theories of 
gravity, in particular general relativity, are ex¬ 
tensively detailed in the monograph by Iwiiil 
(1993). In particular , for f(R) metric theories, 
ICapozziello fc Stabilel (120091 ) have developed a 
perturbation an alysis applied to lenses and clus¬ 
ters of galaxies ( Capozziello et al.1l2009f ). 


In this article, Einstein’s summation convention 
over repeated indices is used. Greek indices take 
values 0,1, 2, 3 and Latin ones 1,2,3. In spherical 
coordinates (x°, x 1 , x 2 , x 3 ) = (ct,r, 0, </?), where 
c is the speed of light, t is the time coordinate 
and r the radial one, with 6 and p the polar and 
azimuthal angles, respectively. The angular dis¬ 
placement dfl 2 := d 6 2 + sin 2 0 dp 2 . We use a 
(+,—,—,—) signature for the metric of the space- 
time. 


Let us consider a fixed point-mass M at the cen¬ 
ter of coordinates generating a gravitational field. 
In this case, the spacetime is static and its spher¬ 
ically symmetric metric g^ is generated by the 
interval 


ds 2 = dx^dx u = goo c 2 dt 2 + #ndr 2 — r 2 dfl 2 , 

(i) 

where g^ is the metric tensor and, due to the 
symmetry of the problem, the unknown functions 
goo and #n are functions of the radial coordinate 
r only. 

The geodesic equations are given by 


d 2 x a dx M dx v 

ds 2 ^ ds ds 


( 2 ) 


where T^ are the Christoffel symbols. In the 
weak field limit when the speed of light c —»• 
00 , ds = c d£, and since the velocity v <C c, 
then each component v l dx°/dt with v l := 
(dr/dt, rdO/dt, rsmOdp/dt). In this case, the 
radial component of the geodesic equations 
for the interval is given by 


In this section we define the relevant proper¬ 
ties of the perturbation theory for applications 


1 d 2 r 
c 2 dt 2 
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where the subscript ( ) ?r := d/d r denotes the 
derivative with respect to the radial coordinate. 

In this limit, a particle bound to a circular orbit 
about the mass M experiences a centrifugal radial 
acceleration given by equation ©, such that: 


u 11 

0 # #oo,r 5 
r 2 


(4) 


for a circular or tangential velocity v. At this 
point, it is important to note that the last equation 
is a general kinematic relation, and does not intro¬ 
duce any particular assumption about the specific 
gravitational theory. In other words, it is com¬ 
pletely independent of the field equations associ¬ 
ated to the structure of spacetime produced by the 
energy-momentum tensor. 


In the weak field lim it of th e th eory, t he metric 
takes the form (see e.g. Landau fc Lifshitz! 1975 b 


#oo — 1 H— 
c z 

922 = ~r 2 , 


- 1 _L 

#11 — —1 H-2~, 

c z 

533 = ~r 2 sin 2 9 , ( 5 ) 


for a Newtonian gravitational potential and an 
extra gravitational potential i/k At the weakest or¬ 
der of the theory, the motion of material particles 
is described by the pote ntial <p, ta kin g int o account 
the fact that 1 / = 0 (e.g. lLandau fe Lifshitzll 19751 ). 
The motion of relativistic massless photon parti¬ 
cles is described by taking into consideration not 
only the second order corrections to the potential 
<fi, but also the sam e order perturbations of the 
potential ip (cf. IWilll 1993 b 


For circular motion about a mass M in the 
weak field limit of the theory, the equations of mo¬ 
tion are obtained when the left-hand side of equa¬ 
tion (j3j) is of order v 2 /c 2 and when the right-hand 
side is of order (p/c 2 . Both are orders 0(c ~ 2 ) of 
the theory, or simply 0(2). When lower or higher 
order corrections of the theory are introduced we 
use the notation 0(n) for n = 0 , 1 , 2 ,... meaning 
O(c 0 ), 0(c -1 ), 0(c ~ 2 ),..., respectively. 


The extended regions of clusters of galaxies 
need a huge amount of dark matter to explain 
the velocity dispersions observed for the stars and 
gas in these systems. At these regions, the veloc¬ 
ity dispersions are typically of order 10 -4 — 10 -3 
times the speed of light. Hence, the “Newtonian” 
physics given by an 0(2) approximation should be 


extended to post-Newtonian 0(4) corrections or, 
equivalently in our model, “post-MONDian” dy¬ 
namics. 

In order to test a gravitational theory through 
different astrophysical observations (e.g. the mo¬ 
tion of material particles, the bending of light- 
massless particles, etc.), the metric tensor g^ v is 
expanded about the flat Minkowski metric for 
corrections h^ u <C r ] flv in the following way: 


9iiv — T h 


(6) 


The metric g^ v is approximated up to second 
perturbation order 0(2) for the time and radial 
components and up to zeroth order for the an¬ 
gular components, in accordance with the spheri¬ 
cal symmetry of the problem. A t this lowest per¬ 
turbation order, Mendoza et ajj ( 2013 ) found the 
time g{^ and radial metric components, for 
the f(x) = X 3 / 2 metric theory of gravity. Our 
notation is such that the superscript (n) denotes 
the order 0(n) at which a particular quantity is 
approximated. These metric values are necessary 
to compare with the astrophysical observations of 
the motion of material particles and that of mass¬ 
l ess p hoton partic les through the bending of light 
(Iwiiillimi Eoofih . in fact, through the observa¬ 
tions of the rotation curves of galaxies and the 
Tully-Fisher relation, and the details of the gravi¬ 
tational lens ing in individua l , grou ps and clusters 
of galaxies, iMendoza et al.l ( 2013 ) fixed the un¬ 
known potentials (p and ip of the theory. 

In this paper, we develop perturbations of the 
relativistic extended model f(x) = X 3 ^ 2 up to the 
fourth order in the time component g$, corre¬ 
sponding to the next order of approxi matio n to de¬ 
scribe the motion of massive particles (jWilll Il993h . 
In this case, the metric components can be written 


as 


#00 —1 + 9oo + goo + 0(6), 

#H = “ 1 + #1? + e) ( 4 )^ /rA 

(0) 2 ^ 

#22 =#22 = ~r , 

#33 =#22 Sin 2 0. 

In other words, the metric is written up to the 
fourth order in the time component and up to the 
second order in the radial one. The contravariant 


4 





















metric components of the previous set of equations 
are given by 


g 00 — 1 “ 9m - 9oo + ^( 6 )> 
5 n = -i -g® 
g 22 =g 22{0) 


0(4), 
1/r 2 , 


(8) 


and is a coupling constant of order one calibrated 
through astrophysical observations. 

Equation m is a particular case of a full 
gravity-field action formulation in which the de¬ 
tails of the mass distributi on appear inside the 
gravitational ac tion (see e.g. lCarranza et ahll2013 : 
Mendoza 20151 for further investigation). In the 


g 33 —g 22 / sin 2 0. 

3. Extended /(x) = y 3 / 2 gravity 
3.1. Field equations 

The f(x) metric theory, proposed bv lBernal et al, 


(1201 lbh . is constructed throu gh the inclusion of 


MOND’s acceleration scale ao ( Milgromlll983b ) as 


a fundamental physical constant, which has been 


present work, we assume the solution for the grav¬ 
itational potential of a point-mass source and 
generalize the result to a mass distribution, for 
applications to spherically symmetric systems, in 
particular clusters of galaxies. 

Now, the matter action takes its ordinary form 

Sm = -1 J Cm x^gd 4 x, (13) 

with £ m the Lagrangian matter density of the sys- 


vance 

3 (see e.g. Bernal et al. 2011a: C 

!arr 

anza et al. A 

2013; 

Mendoza et al. 201 it Mendozal 

2012t Hernandez et ab? 

2010, 

2012: Hernandez & Jimene^ 2012[ 

Mendoza et al. 

2013; 

Mendoza & Olmo* 2014; iMendoza 2015^. 




The correct dimensional Hilbert’s action in 
the metric approach, for a point-mass source 
generating the grav itational field, is given by 
( Bernal et al. 1201 lah 


/'(x)x^ - 3. f(x)9- L 2 m (V„V„ - g^A) f'( x ) 


8t tGL\ 


Mr 


L Hi”) 


(14) 


ft = 


I67 tGL 2 m 


J fix) V^gd 4 x, 


(9) 


for any arbitrary dimensionless function /(x) of 
the dimensionless Ricci scalar 


X l 2 m r , 

where R is the standard Ricci scalar and 

Lm :=((r g l M ) 1/2 , 

is a length scale where 


( 10 ) 


( 11 ) 


where the prime denotes the derivative with re¬ 
spect to the argument, the Laplace-Beltrami op¬ 
erator is A := and the energy-momentum 

tensor is defined through the standard relation 
SS m = —(1/2 c)T a p5g al3 . Also, in equation ([T4]h 
the dimensionless Ricci tensor is defined as 


X/jlv L 2 M R l 




(15) 


where R^ u is the standard Ricci tensor. The trace 


GM 


TV : = 


Im •= 




a 0 


( 12 ) 


with r g the gravitational radius of the sys¬ 
tem, Im the mass-le ngth scale of the system 
( Mendoza et al.1 lioil h a 0 := 1.2 x 10 -10 m/s 2 
th e Milg ro m’s accelerati on constant (see e.g. 
Famaev & McGaugh 2_012, and references therein) 


of equations (fl4l) is given by 


/ / (x)x~2/(x) + 3L| f Af(x) = 8?r / 4 LM r, (16) 
where T := T a a . 

Bernal et al.1 ( 2011bl ) and lMendoza et all (2013) 


have shown that the function /(x) must satisfy the 
following limits: 


r, x _ J X, when x > 1, ( . 

f\X) ^ 3/2 i -1 Gw 

lx', when x 1, 
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in order to recover Einstein’s general relativity 
in the limit y 1 and a relativistic version of 
MOND in the regime % < 1. In the latter case, 
the first two terms on the left-hand side of the 
trace m. are much smaller than the third one 


( Bernal et al. 2011 b), and so 


/ , (x)x-2/(x)«3L 2 m A/ , ( X ), (18) 

at all orders of approximation. All these facts 
mean that the trace m is given by: 

3L 2 M Af\ X ) = (19) 

For the field produced by a point mass M, the 
right-hand side of equation m is null far from 
the source and so, the last relation in vacuum can 
be rewritten as: 


A f{x) = 0- (20) 

Now, for simplicity, we assume a power-law form 
for the function /(%): 

fix) = x b , ( 21 ) 

for a real power b. In this case, relation flUD is 
equivalent to 


A f\R) = 0, (22) 

at all orders of approximation for a power-law 
function of the Ricci scalar 


needed to yield the correct rela tivistic extension 
of MOND. Mendoza et al.l ( 2013 ) showed that this 
condition is crucial to describe the details observed 
for gravitational lensing for individuals, groups 
and clusters of ga laxies, an d differs from the re¬ 
sults obtained bv ICapozziello et al.l ( 2007 ). for a 
standard f(R) power-law description in vacuum. 
As discussed in Appendix A of Mendoza et al. 


( 2013 ), such discrepancy occurs from the sign con¬ 
vention used in the definition of the Riemann ten¬ 
sor, giving two different choices of signature that 
effectively bifurcate on the solution space, a prop¬ 
erty which does not appear in Einstein’s general 
relativity. This is due to higher order deriva¬ 
tives with respect to the metric tensor that appear 
on metric theories of gravity (cf. equatio ns (fl4l) 
and (Il6| ) ). Following the results bv lMendoza et al 


(2013b we use the same definition and branch of 


solutions which recover the correct weak field limit 
of the theory to explain the rotation curves of spi¬ 
ral galaxies, based on the Tully-Fisher relation, 
and the gravitational lensing observed in groups 
and clusters of galaxies, at the outer regions of 
these systems. 

The standard perturbation analysis for f(R) 
metric theories restricted by the constraint equa¬ 
tion 4221 ) is developed for a power-law descrip¬ 
tion of gravity in the weak field limit and in the 
first MOND-li ke relativistic correctio n (cf. equa¬ 
tion (fT7]) h bv iMendoza et al.l ( 2013 ). The stan¬ 
dard field eq uations (fl4l) can be written as (see 
e.g. ICapozziello fe Faraoni 2011) 


f(R) = R b . (23) 

Substitution of the power-law function m in 
the null variations of the gravitational field’s ac¬ 
tion (|9|) in vacuum leads to 

SSt = “ 1 J R* 1 V~ 9 dfx = 0 , (24) 

and so 

sjR b ^x = 0. (25) 

From the last relation, the same field equations 
as the null variation of the action for a standard 
power-law metric f(R ) theory (|23|) in vacuum are 
obtained, but with the important restriction (f22l) 


f\R)R^ - \f{R)g^ + H^ = o, (26) 

where the fourth-order terms are grouped into the 
term 

Rfiv := - (V M V v - A) f\R) . (27) 

The trace of equation 429 is given by 

f'(R)R-2f(R) + H = 0, (28) 

with 

n := U^ v = 3A f'(R) . (29) 
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For the case of the static spherically symmet¬ 


ric spacetime it follows that (iMendoza et al 


201 


n 


Hv 


f r ^R^r 9nv 


11 

,r 


+ 9 11 ( ln V-9) R,r + g^Rrr 

-/"'[R.R.-g^R^Y 


} 


(30) 


and 


n = 3/" 


(fl U r + 9 11 (In y/-9),r) R ,r + l 


5 11 -R, 1 


+ 3 f"9 X1 R 2 r . 


(31) 


The general field equations are of fourth order 
in the derivatives of the metric tensor. In dealing 
with the algebraic manipulations for the perturba¬ 
tions to our f(R) theory of gravity, T. Bernal, S. 
Mendoza and L.A. Torres have developed a code 
in the Computer Algebra System (CAS) Maxima, 
the MEXICAS (Metric EXtended-gravity Incor¬ 
porated through a Computer Algebraic System) 
code (licensed with a GNU Public License Ver¬ 
si on 3). The code is d escribed in Appendix B 
of Mendoza et al. (12013 ) and can be downloaded 
from: http://www.mendozza.org/sergio/mexicas, 
In a more general description of the grav¬ 
itational field, the mathematical form of the 
field’s action (j9j) includes the Schwarzschild mass 
(through Lm) into the integral. This means a 
modification to the standard Lagrangian descrip¬ 
tion of the gravitational field since the matter 
content is generally assumed to appear only in the 


matte r a ction dm Howev er, accord ing to Sobo uti 


( 2007 b iRosas-Guevara! ( 20061 b IMendoza et al 


( 2013 ). modifications at the very fundamental 
level of Hilbert’s action could be expected when 
discussing extensions of gravity. 


3.2. Weakest field limit and post-MONDian 
correction 

We present in this subsection fourth order per¬ 
turbation results for the metric coefficients and the 
Ricci scalar. 


Ricci’s scalar can be written as follows: 

R = i? (2) + R {4) + 0(6), (32) 

which has a non-null second and higher perturba¬ 
tion orders. This scalar is a function of the met¬ 
ric components and their derivatives with respect 
to the coordinates up to the second order. The 
fact that R^ = 0, is consistent with the flatness 
of spacetime assumption at the lowest zeroth per¬ 
turbation order. The second order component of 
Ricci’s scalar R^ from the metric components ([7]) 
is given by 



( 2 ) 

a (2) + ^- 
9ll,r + r 


(2) 

#00, rr 


2 ( 2 ) 
r #00 ,r > 


(33) 


and the fourth order component is: 


r { a) 


$ 


2 #ll!r + #00 ]r 


( 2 ) 

#11 


( 2 ) 

' #00,r 




(2) A 

-#00,rj 


( 2 )' 

#00 


r ( 2 ) 

o #00,rr 


( 2 ) ( 2 ) 
#00 - #11 


Z (4) (4) 

^#00,r #00,rr ‘ 

(34) 


At the lowest perturbation order 0(26 — 2), the 
trace of the field equations for a theory given by 
equation ([23]) . can be written as (IMendoza et al 


2013D: 


^(26-2) = 3A/ , ( 26 - 2 )( j R) = 0. 


(35) 


Note that this is the only independent equation 
at this perturbation order. Substitution of ex¬ 
pressions (0, ([23]) and ([32]) in the previous equa¬ 
tion leads to a differential equation for Ricci’s 
scal ar at order 0(2), wh ich has a solution given 
by ( Mendoza et al. 2013 ) 


i? (2) (r) 



i/(6— i) 


(36) 


where A and B are constants of integration. Far 
away from the central mass spacetime is flat, and 
so, Ricci’s scalar must vanish at large distances 
from the origin, i.e. the constant B = 0. The 
case 6 = 3/2 yields a MONDian-like behavior for 
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the gravitational fie l d in the limit r I m r g 
( Bernal et al • 2011 bMMendoza et alJll2Q13h and so, 
substituting 6 = 3/2 in relation ([36]) and 23 = 0 
yields: 


R {2 \r) = 5, 


(37) 


where R := A 2 / 4:. 

At the next perturbation order 0(26), the met¬ 
ric components g^, an d Ricci’s scalar 

R can be obtained. The field equations (]26l) at 
this order are given by (Mendoza et al. l2013f) 


bR^-'Rff - = 0, (38) 

where 

= - (V„ V„ - A) /'( 2b ) (i?). (39) 


and using equation dS3, the following differen- 

( 2 ) - 

tial equation for g\^ is obtained ([Mendoza et al 


(2) , ( 2 ) , R _ n 

r 9ll, r + 511 + ~2 + 77 “ 


with analytic solution 


(44) 


(2)^ + 


flirw = 


3 ’ 


(45) 


where /c 2 is a constant of integration. 

Now, to derive the first relativistic correction of 
the metric theory (f2Tj) with b = 3/2, i.e. to obtain 
g^ and R^ having in mind further applications 
to the outer regions of clusters of galaxies, we use 
another independent field equation, namely the 22 
component of equation ([38]) . In fact, substitution 
of relations © and ([8]) in the 22 component of 
equation ([38]) yields 


From equation ([22]) it follows that A/'( 26 ) = 0 
and the last equation simplifies greatly. Using re¬ 
lations ® and (Et in the 00 componen t of equa¬ 
tion (]38l) leads to ( Mendoza et al.ll2Q13h 


b(b - 1 )rR^ b ~ 2 'rW + g[^R (2) + (b - 2 )rW~ 1 rWrW 


- bR.Wh-'-R^ - — R (2)b = 0, 


(46) 


fo J RC 2 )b-l jR < 2 )_l jR ( 2 ) 6 _ | _l b ( b _i)^) pjR ( 2 ) 6 - 2 j i C 2 ) = 0> 

(40) 

where the 00 component of the Ricci tensor at 
0(2) is 


rn ( 2 ) i 9 

(2) _ ' VQQ,rr ^ z f/QQ,r 


; — — 

^00 


2r 


(41) 


Substituting this last expression, 6 = 3/2 and re¬ 
sult ([37]) into equation ([40]) , the following differen¬ 
tial e quation for is obtained ( Mendoza et al 
2013h : 


where the 22 component of the Ricci tensor at 
0(2) is given by 




( 2 ) 

9oO,r 


011,r 


(47) 


Using the fact that 6 = 3/2, and substituting the 
last equation together with relations ([37]) and ([45]) 
into equation d®, we obtain the following differ¬ 
ential equation for the Ricci scalar at 0(4): 


r 4 R < 4) + r 3 i? (4) + R 2 r - 3 k 2 R = 0, (48) 

which has the following exact solution: 


2 ( 2 ) 
r 000 rr 


(2) 2 R 

+ 3r 5oO,r + “V = 


which has a solution given by 


(42) 


5ooW = —+^» («) 

where k\ and r s are constants of integration. 
By substitution of this result into equation (133]) 


i? (4) (r) 


R 2 3 k 2 R 4fc 3 
r 2 2 r 3 r 


(49) 


where k% is a constant of integration. Now, the ex¬ 
pression for the Ricci scalar R from the metric 
components © is given by expression ([Ml)- Sub¬ 
stituting such equation in ([49]) , together with rela¬ 
tions d®, d® and d®, we obtain the following 
differential equation for g^: 
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-9 g 


(4) 

00 , rr 


18 

r 


(#00 ]r + 




23 

~2~ 


r A \ r J r 4 

45/ci/c 2 54fc? 

+ r 5 + ^6“ = °> 

(50) 



3.3. Generalization to extended systems 

To compare the correction g qq with the obser¬ 
vations of clusters of galaxies, let us take the radial 
component m of the geodesic equations ([2]) in the 
weak field limit of the theory. In this limit, the ro¬ 
tation curve for test particles bound to a circular 
orbit about a mass M with circular velocity is v(r) 
given by equation Such equation, up to the 
fourth order of approximation, is given by 


with the following exact solution: 


( 4 ) _^ k 2 

%o ~ 18 m 




7 2/^4 

2 fa r-b fa 

r 


In ( — 


fa R 
3 r 2 


fa 1 5 fa fa fa 2 

j— t —|-I—1 _i—_ 

7*2 6 2 r 4 5 


hfaR 

3r 



(51) 


where k 4 and fa are constants of integration. 

To fix t he free param e ters in relations (|43]) 
and IMendoza et al.l (I2013I) compared the 

metric coefficients with observations of rotation 
curves of spiral galaxies and the Tully-Fisher re¬ 
lation and with gravitational lensing results of in¬ 
dividual, groups and clusters of galaxies. Their 
results are summarized in Table [lj 


Using these solutions of the metric coefficients 
at order 0(2), the metric coefficient g$ and 
Ricci’s scalar R ^ reduce to 



(52) 

for v z /Im = (GMao) 1 / 2 /c 2 . 

To fix the constants of integration fa and fa (fa 
vanishes upon derivation of with respect to r 
- cf. equation ([53]) h we adjusted the observational 
data of 12 clusters of galaxies, as described in the 
next subsection. 


1 d 2 r 
c 2 d t 2 


( 2 ) 

# 00,r 


( 2 ) ( 2 ) 

' #11 #00,r 


(4) 

" #00,r 


(53) 


Substitution of the 0(2 ) perturbation values 
of the metric coefficients from Table [T| and equa¬ 
tion (l5T1) in ([53]) , results in the following expression 
for the acceleration of a test mass particle in the 
gravitational field generated by the point-mass M: 


V) = 


(GMa 0 ) 1/2 


r 

23GMa 0 


2 GMa 0 


In 


+ c 4 k 3 - 


c 4 &4 


(54) 


where a c := \a c \. The first term on the right-hand 
side of last equation correspon ds to the u deep- 


MONDian” acceleration (see e.g. lFamaev fe McGaugfj 
20120 . The remaining 0(2) terms are the first 
relativistic correction to the gravitational acceler¬ 
ation. 

In order to apply these results to clusters of 
galaxies, it is necessary to generalize the gravita¬ 
tional acceleration ([54]) to a spherical mass distri¬ 
bution M(r). The first term of such equ atio n can 
be eas ily generalized according to IMendoza et al. 


( 201 lh : In this case, the deep-MONDian accel¬ 
eration can be written as f(x) = a/ap = x, 
for x := Im/t. As discussed bv IMendoza et all 
(|2011 ) this function, and in general any analyti¬ 
cal function which depends only on the parameter 
x , guarantees Newton’s theorems. In other words, 
the gravitational acceleration exerted by the outer 
shells at position r cancels out and depends only 
on the mass M(r) interior to r: 


M(r) = 47r / p(r)r 2 dr, 

J 0 


(55) 
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Table 1 


Metric coefficient 

( 2 ) 

%0 

o (2) 

9 11 


-r^in (f 1 ) 

l M \ r s J 

_2rg 

Observations 


(Tully-Fisher) 

(Lensing) 

Theory 

-#ln(*) + ^ 

k l \ k 2 R 

7F + - - ~3 

fix) = x 3/2 


R = 6v s /Im , ki = 0 

R = 6v s /Im, &2 = 0 


Note. —The table reproduced from lMendoza et al.l i2013h shows the 
results obtained for the metric components g^ and g[ 2 ± for a static 
spherically symmetric spacetime. The metric coefficients are empiri¬ 
cally obtained from astronomical observations in the scales of galaxies 
(Tully-Fisher relation) and lensing in the outer regions of individual, 
groups and clusters of galaxies, and compared to the ones predicted 
by the metric /(%) — X 3//2 theory of gravity. Any proposed metric 
of a theory of modified gravity must converge to the observational 
values presented in the table. As shown by the authors, the theory 
fix) = X 3/2 is i n perfect agreement with the observational metric 
components. Since the metric components up to the order 0(2) deter¬ 
mine the “MONDian gravitational potential” of the system, the length 
r s is undetermined. However, as explained in section^ its value is nec¬ 
essary to describe the dynamics up to perturbation order G( 4) of the 
theory and so it will be necessary to fix its value with observational 
data. 


where p is the mass density of the system. Thus, 
the first term of the gravitational acceleration ([54]) 
due to a mass distribution can be written as 


is composed of many infinitesimal mass elements 
dm, each one contributing with a point-like grav¬ 
itational potential (EZD, such that: 


r 

For the second order 0(2 ) terms on the right- 
hand side of acceleration (EH), let us take the grav¬ 
itational potential generated by a point-mass par¬ 
ticle M for the model f(x) = X 3 • Integrating 
the acceleration ([54]) with respect to the radius r, 
according to the spherical symmetry assumption, 
we obtain 


m\ 



dr, 


+ 


1 


= (GMa 0 ) 1/2 ln(^j 


|cMao In 



23-In 



+ c A ksr + 



(57) 


This point-mass gravitational potential can be 
generalized considering that the extended system 


M(r) = [ dm = [ p(r')dV', (58) 

Jv Jv 

where the volume element is dV' = r' 2 sin 6' d(p' d O' dr', 
integrated over the volume V. Notice that the 
mass of the system M appears elevated to the 
one half power in the first term of equation dsn, 
and that the potential depends linearly with the 
mass in the second and third terms. Since the 
constants k% and k 4 are proportional to 1/c 4 in 
order to have second order terms in the acceler¬ 
ation, let us assume that they are proportional 
to GMao/c 4 . Thus, the 0(2) point-gravitational 
potential can be written as 


. (j)^ 2 \r) 


GMa 0 



23-In 



+ Ar — 
(59) 


where A and B are constants. 

From equation (|531i . the generalized gravita- 
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tional potential <f>( 2 )(r) in spherical symmetry is 
the convolution 


J /(r — r')p(r')r' 2 sin 6'dp'd6'dr', (60) 


of the function 


( Ga, 

/ r-r = — 

c z 




23-In 


B 


(61) 


+A\r — r'| - 


with the differential dm, given by equation ^ 
(see e.g. IVladimirov 2002 ). for / and p locally in¬ 
tegrate functions for r > 0. Due to the spheri¬ 
cal symmetry of the problem, the integration can 
be done in one direction, for example the z axis, 
where the polar angle 6 = 0 and | r — r' | = 
\Jr 2 + r' 2 — 2 rr' cos O' . Thus, the 0(2) general¬ 
ized gravitational potential for a mass distribution 
can be written as 


-w-Wf {'• 


r — r 


23-In 


r — r 


+A\r — r'\ 


) r' 2 sin#' dp' d 6' dr', 
(62) 


integrated over the whole volume. If the den¬ 
sity distribution is known, the extended poten¬ 
tial ([62]) can be numerically integrated to obtain 
the gravitational acceleration, from 0 < r < r' and 
r' < r < 7Z, where 7 Z is the radius of the spherical 
configuration. It is a well-known result that the 
matter outside the spherical shell of radius r does 
not contribute to the potential for the fourth term 
on the r ight-hand side of equation (l62j) (see e.g. 
Vl adimi rov 120020 . For the other terms, the inte¬ 
gration is done for the interior and exterior shells 
of mass with respect to the radius r, giving as re¬ 
sult the following expression: 


$( 2) (r) 


27rGao 




r — r 


{r + rVln l±IL 


In ■ 


1 r + r' 

- m- 

2 r s 


1 | r-r' 

- m- 

2 r. s 


24 


(r + r') 3 - |r - 


j 13 


3 L 

12Ga 0 rU 


| p(r')r' dr' 


/ 47rp(r') 

Jo 


r' 2 dr', 


(63) 


where the last term is constant. Performing the 
derivation of the potential ([62]) with respect to r 
and simplifying some terms, the generalized grav¬ 
itational acceleration can be written as 


_ _ [GM{r)a 0 ] 1/2 , d$( 2 )(r) 

a C\ r ) ~ I jj 5 

r dr 

_ [GM(r)ao] 1/2 GMipapB d$j 2) (r) 

r c 2 r 2 dr ’ 

(64) 

which can be obtained for a given p(r'). Notice 
that the parameters A and r s appear only on the 
last term of equation ([64]). 

4. Cluster mass profiles 

4.1. Fit with observations of clusters of 
galaxies 

To apply the results of the last subsection to 
the spherically symmetric X-ray clust ers of galax¬ 
ies reported by IVikhlinin et al.l ( 2006 ), notice that 
there are two observables: the ionized gas profile 
p g (r) and the temperature profile T(r). Under the 
hypothesis of hydrostatic equilibrium, the hydro- 
dynamic equation can be derived from the colli¬ 
sionless isotropic Boltzmann equation for spheri¬ 
cally symmetric systems in the weak f ield limit of 
approximation (see e.g. Binnev fc Tremainell2008l) : 


(65) 

where <f> is the gravitational potential and cr r , cfq 
and dp are the mass-weighted velocity dispersions 
in the radial and tangential directions respectively. 
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For an isotropic system with rotational symmetry 
there is no preferred transverse direction, and so 
(tq = dcp. For an isotropic distribution of the ve¬ 
locities, we also have a r = cfq. 

The radial velocity dispersion can be related 
to the pressure profile P(r), the gas mass density 
p g (r) and the temperature profile T(r) by means 
of the ideal gas law to obtain: 


.2 = = fc B T(r) 

r Pg(r) P m p 


( 66 ) 


where k b is the Boltzmann constant, fi w 0.609 
is the mean molecular mass per particle and m p 
is the mass of the proton. Direct substitution 
of equation ([66]) into ([65]) yields the gravitational 
equilibrium relation: 


d $(r) 

k B T(r) 

dlnpg(r) | din T(r) 

dr 

lim p r 

d In r d In r 


(67) 

The right-hand side of the previous equation is 
determined by observational data, while the left- 
hand side should be consistent through a given 
gravitational acceleration and distribution of mat¬ 
ter. In standard Newtonian gravity, the total mass 
of the cluster is given by the mass of the gas and 
the mass of the galaxies inside it, with the ad¬ 
dition of an unknown dark matter component to 
avoid a discrepancy of one order of magnitude on 
both sides of the last equation. 

For the particular Newtonian gravity case, the 
“dynamical” mass Md yn of the system is deter¬ 
mined by the Newtonian acceleration a n: 


M — r ^ N - feBT ( r ) d In p g (r) d In T(r) ~ 

dyn * G fim v G T dlnr dlnr 

( 68 ) 

In the /(x) = x 3 ^ 2 model, the acceleration will 
be given by equation ([6l]h thus we can define the 
“theoretical” mass M t h as: 


M t h : = 


r 2 a c (r) 

G 


~M{r)a 0 1 1/2 M(r)a 0 B r 2 d<f>¥\r) 

G J r+ t? + G dr 

(69) 


where the mass M(r) is given only by the baryonic 
mass in the system, i.e.: 


M(r) = M gas (r) + M stars (r). (70) 

In the previous equation, M gas is the mass of the 
gas and M sta rs is the stellar mass in the cluster. 

In order to reproduce the observations, the 
theoretical mass obtained from our modification 
to the gravitational acceleration must be equal 
to the dynamical one coming from the observa¬ 
tions, without the inclusion of dark matter, i.e. 
M t h = Mdyn- This provides an observational pro¬ 
cedure to fit the three free parameters r s , A and 
B of our model. 


4.2. Chandra clusters sample 


Chandra X-ray clusters of galaxies sample can 
be found in the article of IVikhlinin et al~ ( 2006f ) 
and include the full analysis of 12 clusters: gas 
density, temperature and total mass profiles. 

The observed keV temperatures in clusters of 
galaxies are interpreted in such a way that the 
gas is fully ionized and the hot plasma is mainly 
emitted by free-free radiation processes. There 
is also line emission by the ionized heavy ele¬ 
ments. The radiation process generated by these 
mechanisms is proport ional to the em is sion m ea- 
sure profile n p n e (r). Vikhlinin et al.l (|2006f ) in¬ 
troduced a modification to the sta ndard /3-model 
of ICavaliere &; Fusco-Femianol ( 1978 b in order to 
reproduce the observed features from the surface 
brightness profiles, the gas density at the centers of 
relaxed clusters and the observed X-ray brightness 
profiles at large radii. A second /3-model compo¬ 
nent (with small core radius) is added,jto inc rease 


accur acy near the cluster centers ([Vikhlinin et al 


2006). With these modifications, the complete ex¬ 
pression for the emission measure profile has 9 free 
parameters. The 12 clusters can be adequately fit¬ 
ted by this model. The best fit values to the emis¬ 
sion measure for th e 12 clusters of galaxie s can be 
found in Table 2 of IVikhlinin et all ( 2006 ). 

To obtain the baryonic density of the gas, 
the primordial abundance of He and the relative 
metallicity Z = 0.2 Zq are taken into account and 
so, 


p g (r) = 1.624ra p 



(71) 
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For the stellar component of the clusters, we 
used the empirical relation between the stellar and 
the total mass (baryonic plus dark m atter ) in the 
Newtonian approximation (Lin et al. 2012): 


M si 


10 12 M© 


= (1.8 zb 0.1) 


-^500 

10 14 M© 


0.71±0.04 


, (72) 


where M500 is the mass of the cluster in the stan¬ 
dard dark matter description. 


For the temperature profile Tfrb lVikhlinin et al.l 
( 2006 ) used a different approach from the poly¬ 
tropic law, i.e. T(r) oc Pg -1 , to model non¬ 
constant cluster temperature profiles at large 
radii. All the projected temperature profiles show 
a broad peak near the centers and decreases at 
larger radii, with a temperature decline toward 
the cluster center, probabl y because of the pres - 
ence of radiative cooling ( Vikhlinin et al. 2006). 
To model the temperature profile in three dimen¬ 
sions, they constructed an analytic function such 
that outside the central cooling region, the tem¬ 
perature profile can be represented as a broken 
power law with a transition region. The 8 best-fit 
parameters for the temperature profiles of the X- 
ray clusters of galaxie s can be found in Table 3 of 
Vikhlinin et al. ( 2006 ). 

The total dynamical masses obtained with the 
use of equation (|68|) from the derived gas den¬ 
sity and temperature profiles for 12 Chandra X- 
ray clusters of galaxies, were kindly provided by 
A. Vikhlinin. 


4.3. Parameters estimation method 

We conceptualized the three parameters cali¬ 
bration, A , B and r s , as an optimization problem 
and propose to resolve it using Genetic Algorithms 
(GAs), which are evolutionary based stochastic 
search algorithms that, in some sense, mimic nat¬ 
ural evolution. In this heuristic search technique, 
points in the search space are considered as in¬ 
dividuals (solution candidates), which as a whole 
form a population. The particular fitness of an in¬ 
dividual is a number, indicating their quality for 
the problem at hand. As in nature, GAs include 
a set of fundamental genetic operations that work 
on the genotype (solution candidate codification): 
m utation , recombination and selection operators 
riMit,chelllll998h . 


These algorithms operate with a population of 
individuals P(t) = x^,...,x^, for the t-th iter¬ 
ation, where the fitness of each X{ individual is 
evaluated according to a set of objective functions 
fj{xi). These objective functions allow to order, 
from best to worst, individuals of the population in 
a continuum of degrees of adaptation. Then, indi¬ 
viduals with higher fitness, recombine their geno¬ 
types to form the gene pool of the next generation, 
in which random mutations are also introduced to 
produce new variability. 


A fundamental advantage of GAs versus tradi¬ 
tional methods is that GAs solve discrete, non con¬ 
vex, discontinuous, and non-smooth problems suc¬ 
cessfully, and thus they have been widely used in 
Ecology, Natural Resources Management, among 
other fields, but no t so much in Astrophysics 
( Lopez-Corona et al.l 20131) . Neverth e less, they 
have been recently used by iNesseri 1 (l201ll) for 
parameters searches in ACDM models with SNIa 
data. 


It is important to note that, as it is well known 
from Taylor series, any (normal) function may be 
well approximated by a polynomial up to certain 
correct order of approximation. Of course, even 
that this is correct from a mathematical point of 
view, it is possible to consider that a polynomial 
approximation is not universal for any physical 
phenomenon. In this line of thought, one may fit 
any data using a model with many free parame¬ 
ters, and even in this approximation we may have 
a great performance in a statistical sense, but it 
could be incorrect from the physical perspective. 


In this sense, an important question to ask is, 
how much better a complex (more parameters) 
model must perform in a fitting process, justifying 
the incorporation of extra parameters? In a more 
straightforward sense, how do we trade off fit with 
simplicity? Such question has been the motivation 
in the recent years for new model selection criteria 
development in statistics, all of which define sim¬ 
plicity in terms of the number o f parameters, o r 
the d imension of a model (see e.g. IForster fc Sober 
1 1994 for a non-technical introduction). These cri¬ 
teria include Akaike’s Information Criterion (AIC) 
( Akaikel Il974 Il985l). the Bay esian Information 
Criterion (BIC) ( Schwarz1ll978) and the Mi nimum 
Description Length (MDL) ( Rissaneri19891 ). They 
fit the parameters of a model a little different be¬ 
tween them, but all of them address the same 


13 
























































problem as a significance test: Which of the esti¬ 
mated “curves” from competing models best rep¬ 
resent reality? ( Forster fe Sober [l994j h 

Akaikel ( 1974 1985 ) has shown that choosing 


the model with the lowest expected information 
loss (i.e., the model that minimizes the expected 
Kullback-Leibler discrepancy) is asymptotically 
equivalent to choosing a model Mj, from a set 
of models j = 1,2that has the lowest AIC 
value, defined by: 


the average from the first best population decile, 
checking the consistency of the second order cor¬ 
rections with respect to the zeroth order term in 
the gravitational acceleration dSU). Finally, since 
we obtained: 


Aaic := AIC* - min {AIC*} < 2, (74) 

for the parameters estimation, then the m odel was 
accep ted as a good one ( Burnham fe Andersonl 


2002 


AIC = -21n(£ i ) + 2V-, (73) 


5. Results and discussion 


where Cj is the maximum likelihood for the can¬ 
didate model and is determined by adjusting the 
Vj free parameters in such a way that they max¬ 
imize the probability that the candidate model 
has generated the observed data. This equation 
shows that AIC rewards descriptive accuracy via 
the maximum likelihood, and penalizes lack of par¬ 
simony according to the number of free parameters 
(note that models with small er AIC values are to 
be preferred). In that sense, lAkaikel ( 1974 Il985 ) 
extended this paradigm by considering a frame¬ 
work in which the model dimension is also un¬ 
known, and must therefore be determined from 
the data. Thus, Akaike proposed a framework 
where both model estimation and selection could 
be simultaneously accomplished. For those rea¬ 
sons, AIC is generally regarded as the first, most 
widely known and used model selection tool. 

Taking as objective function the AIC informa¬ 
tion index, we performed a Genetic Algorithm 
analysis using a modified version of the ISastrvI 
(2007) code in C++. The GA we used evalu¬ 
ates numerically equation ([69]) in order to com¬ 
pare the numerical results from the_model with 
the observational cluster data by Ivikhlinin et al 


( 200 6). as explained in subsection 14.21 All pa¬ 
rameters were searched in a broad range from 
— lx 10 4 to 1 x 10 10 , generating populations 
of 1,000 possible solutions over a maximum 
of 500,000 generation search processes. We 
selected standard genetic algorithms: tourna¬ 


ment selectio n with replacement (jGoldberg et al 
19891: Sastrv fe Goldberg] 200 ll). simulat e d bi¬ 
nary crossover (SB X) ([Deb fe Agrawall Il995 


Deb fe Kumarl 1 19951) and polynomial mutation 
( Deb fe Agrawall ' 19951 : iDeb fe Kumar 19951: IPebl 


200 1). The parameters were estimated taking 


The results for the best fits as explained in the 
last subsection are summarized in Table O Fig¬ 
ures 1 and 2 show the best fits compared to the 
total d ynamical mass obtained by Ivikhlinin et al. 
( 2006 ). as mentioned in section 0J 

From the best-fit analysis, we see that our 
model is capable to account for the total dynami¬ 
cal mass of the 12 clusters of galaxies, except at the 
very inner regions for some of them, a persistent 
behavior more accentuated for A907 and A1991. 
Furthermore, as can be seen in table [2j there are 2 
clusters, A1413 and A2029, for which the parame¬ 
ter r s estimated is very far from the mean value of 
the other clusters. This parameter appears in the 
gravitational acceleration ([64]) only through the 
integral ([63]) . Comparing its contribution to the 
acceleration with respect to the other two terms 
in equation we found that the dominant sec¬ 
ond order term is the one with the parameter F>, 
and the contribution of the derivative of the in¬ 
tegral ([63]) is very small (because r s appears in¬ 
side a logarithm and the particular combination 
of the functions in such equation). Moreover, the 
Aaic < 2 obtained for all the clusters indicates 
that, in general, the fitted model is very good. 

From the figures, we see that our “MOND-like” 
relativistic correction model is better at the outer 
regions of these systems, exactly where dark mat¬ 
ter is introduced in the standard Newtonian grav¬ 
ity scenario. In this sense, our model is better than 
standard MOND, which needs huge amounts of ex¬ 
tra matter to fit the observations in these systems. 
Also, the second order perturbation analysis of the 
metric theory CD with b = 3/2, was capable to 
account for the observations of the rotation curves 
of spiral galaxies and the Tully-Fisher relation, 
and the gravitational lensing in individual, groups 
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Fig. 1.— Dynamical mass vs. radius for the first 6 clusters of galaxies, with the obtained para meters s umma- 
rized in Table 1. The points with uncertainty bars are the total mass obtained by the fitting of lVikhlinin et al~ 
( 2006 b The solid line is the best fit obtained with our model. 
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Table 2 


Cluster t min (kpc) 

T max (kpc) 

M 6 (1O 13 M 0 ) 

M th (10 14 M o ) 

M dyn (lO 14 M 0 

) A(kpc 1 ) 

B(10 s kpc) 

r s (10 8 kpc) 

A133 

92.10 

1005.81 

3.193 

3.269 

3.359 

- 96.563 

2.0727 

3.42 

A262 

62.33 

648.36 

1.141 

0.825 

0.8645 

- 358.82 

1.7825 

2.40 

A383 

51.28 

957.92 

4.406 

2.966 

3.17 

- 151.99 

1.9714 

2.96 

A478 

62.33 

1347.89 

10.501 

7.665 

8.18 

- 61.436 

1.9271 

3.91 

A907 

62.33 

1108.91 

6.530 

4.499 

4.872 

- 101.05 

1.9024 

3.56 

A1413 

40.18 

1347.89 

9.606 

7.915 

8.155 

- 51.150 

1.9423 

71293 

A1795 

92.10 

1222.57 

6.980 

6.071 

6.159 

- 61.757 

1.9255 

3.79 

A1991 

40.18 

750.55 

1.582 

1.198 

1.324 

- 340.56 

2.3703 

2.49 

A2029 

31.48 

1347.89 

10.985 

7.872 

8.384 

- 75.339 

2.1837 

440.9 

A2390 

92.10 

1415.28 

16.621 

11.151 

11.21 

- 16.935 

1.1517 

4.55 

MKW4 

72.16 

648.36 

0.676 

0.805 

0.8338 

- 367.55 

2.4413 

2.28 

RXJ1159+5531 

72.16 

680.77 

0.753 

1.105 

1.119 

- 171.94 

2.3460 

2.39 






Mean value 

- 154.59 

2.0014 

5980 






< SD > 

1.6462 

0.00868 

937.0 

Note. —Parameters 

estimation for the 12 

clusters of galaxies studied. From 

left to right, the 

columns represent 

the name 

of the cluster, the 


minimal r m i n and maximal r max rad ii for the inte gr atio n, the total baryonic mass Mb, the total theoretical mass M t h derived from our model, the 
total dynamical mass Md yn from fvlkhlinin et al.l (2006), the best-fit parameters A, B and r s , respectively. Also, at the bottom of the table, we 
show the best-fit parameters obtained from the 12 clusters of galaxies data taken as a set of independent objective functions together, with their 
obtained Mean Standard Deviations < SD >. 







Fig. 2.— The same as Figure 1 for the last 6 clusters of galaxies. 
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and clusters of galaxies ([Mendoza et al.ll2013l) . In 
this work, we kept fixed these parameters at per¬ 
turbation order 0(2) to obtain the 0(4) of the 
model, with the additional result that it is possi¬ 
ble to fit the dynamical mass of clusters of galaxies 
without the need of extra dark matter. 

Up to now it has generally been thought that 
a MOND-like extended theory of gravity was not 
able to explain the dynamics of clusters of galax¬ 
ies without the necessary introduction of some 
sort of unknown dark matter component. Our 
aim has been to show that in order to account 
for this dynamical description without the inclu¬ 
sion of dark matter, it is necessary to introduce 
relativistic corrections in the proposed extended 
theory. To do so, we have chosen the particu¬ 
lar fix) = y 3 / 2 MO ND-like metric extension of 
iBernal et al. ( 2011bh . which has also shown to be 
in good agreement with gravitational lensing of in¬ 
dividual, groups and clusters of galaxies and with 
the dynamics of the universe providing an accel¬ 
erated expansion without the introduction of any 
dark matter^and/or energy entities (see e.g. the 
review of iMendozal 12015 . and references therein). 
A similar analogy occurred when studying the or¬ 
bit of Mercury about a century ago. Its motions 
are mostly understood with Newton’s theory of 
gravity. However it was necessary to add rela¬ 
tivistic corrections to the underlying gravitational 
theory to account for the precession of its orbit. 
Mercury orbits at a velocity ^ 50km/s, implying 
a Lorentz factor of ^ 10 -4 and already relativis¬ 
tic corrections are required. Typical velocities of 
clusters of galaxies are ~ 10 3 km/s with a Lorentz 
factor ^ 10 -3 . This means that the dynamics of 
clusters of galaxies are about one order of mag¬ 
nitude more relativistic than the orbital velocity 
of Mercury and so, if the latter required relativis¬ 
tic corrections, then the necessity to describe the 
dynamics of clusters of galaxies with relativistic 
corrections are even more important. 
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